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We study the metastability of the superheated Meissner state in type II superconductors with k ^ 
1 beyond Ginzburg-Landau theory, which is applicable only in the vicinity of the critical temperature. 
Within Eilenberger's semiclassical approximation, we use the local electrodynamic response of the 
superconductor to derive a generalized thermodynamic potential valid at any temperature. The 
stability analysis of this functional yields the temperature dependence of the superheating field. 
Finally, we comment on the implications of our results for superconducting cavities in particle 
accelerators. 
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I. INTRODUCTION 

The Meissner effect - the expulsion of a weak mag- 
netic field from a bulk superconductor - is one of the 
hallmark of superconductivity. As the field is increased, 
the response of a superconductor depends on the value 
of the Ginzburg-Landau (GL) parameter n = A/^, where 
A and ^ are the magnetic field penetration depth and the 
superconducting coherence length, respectively. Type I 
superconductors {k < l/\/2) usually turn normal at the 
thermodynamic critical field He, but superconductivity 
can be maintained, as a mctastable state, up to the super- 
heating field Hsh > He- As for type II superconductors 
{k > l/\/2), above the first critical field Hd their sta- 
ble state is characterized by the presence of vortices, and 
superconductivity persists up to the second critical field 
Hc2- However, since the work of Bean and Livingston^ it 
is known that an energy barrier at the surface impedes 
the penetration of vortices into the bulk, making it pos- 
sible for the Meissner state to exist as a metastable state 
up to a superheating field Hsh > Hci- Therefore Hsh 
is a characteristic property of both type I and II super- 
conductors. In this paper we consider the temperature 
dependence of the superheating field in strong Type II 
superconductors with k 3> 1 - i.e., in the London (or 
local) limit. 

Over the years the issue of the stability of the super- 
heated Meissner state has received much attention. The 
simplest system in which this problem can be studied is a 
clean superconductor occupying a half space with a mag- 
netic field applied parallel to the surface. For this system 
in the strong type II limit, and assuming that the insta- 
bility is due to fluctuations in the direction perpendicu- 
lar to the surface (i.e., one-dimensional fluctuations), de 
Gennes^ calculated the superheating field Hsh = He near 
the critical temperature Tc- If the instability signals the 
penetration of vortices, however, the relevant fluctuations 
can be expected to vary along two dimensions while pre- 
serving translational invariance along the field direction. 
Galaikc^ showed that this is indeed the case and near Tc 



the actual superheating field is smaller than that found 
by de Gennes, Hsh — 0.745i/c- More details about the 
critical fluctuations were presented by Kramer,-^ espe- 
cially in relation to the problem of vortex nucleation. The 
question of metastability has also attracted the interest 
of the mathematical community)^ and a detailed study 
of the instability due to one-dimensional fluctuations in 
type II superconductors was presented in Ref. i7|. A simi- 
lar analysis was performed for type I superconductors in 
Ref. d, in which the results of earlier numerical^iiS and 
analytical^^ investigations are confirmed and extended. 

It is interesting to note that all the previous calcula- 
tions of the superheating field were performed within the 
GL theory (with the exception of Ref. Ji in which the zero- 
temperature limit is also considered). This approach, 
however, is justified only near the critical temperature 
(i.e., if Tc — 2^ ^ Tc), and a quantitative evaluation of 
the superheating field at low temperatures requires the 
use of the microscopic theory of superconductivity. Un- 
derstanding the temperature dependence of the super- 
heating field is important for practical applications: the 
maximum accelerating field of superconducting cavities 
used in particle accelerators is limited by the superheat- 
ing fieldi^ and the optimal operational temperature lies 
wefl below T^}^ 

In this work we consider, within the semiclassical ap- 
proach of Eilenberger,™ a clean type II superconductor 
occupying the half space x > in the presence of an ex- 
ternal magnetic field Ha parallel to the surface. We de- 
rive an expression for the thermodynamic potential valid 
at any temperature for k 3> 1, which enables us to cal- 
culate the temperature dependence of the superheating 
field. As a result we find that in the limit k +oo the 
ratio Hsh/ He between superheating and thermodynamic 
critical fields is a non-monotonic function of temperature 
which has a maximum at T ~ O.OGTc- 

This paper is organized as follows: in Sec. [Til we briefiy 
review the semiclassical theory of superconductivity and 
introduce our notation, while the derivation of the ther- 
modynamic potential is presented in Sec. IIIII The con- 
ditions for the (meta)stability of the Meissner state are 
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discussed in Sec. IIVI and the superheating field is calcu- 
lated in Sec.|Vl In Sec. IVII we discuss the implications of 
our results for accelerator applications and a brief sum- 
mary is given in Sec. IVIII 



II. SEMICLASSICAL THEORY 

In the semiclassical approach to superconductivity the 
superconducting system is described by the set of equa- 
tions, the so-called Eilenberger equations, which are valid 
at any temperature^^ under the assumption that the 
Fermi wave length Af is the smallest length scale char- 
acterizing the system. In practice this means that the 
semiclassical approximation usually applies to low-Tc su- 
perconductors, in which the zero-temperature coherence 
length ^0 ^ Ap. This condition also implies the applica- 
bility at magnetic fields as high as i?c2- This technique 
is widely used to study the properties of hybrid super- 
conducting devices; see, e.g., Ref. [161. From now on, we 
will use units such that Boltzmann constant fc^ = 1 and 
the Planck constant h ^ 1. 

The Eilenberger equations are equations for the 
anomalous Green's functions /(wn, n, r) and /(w„, n, r), 
which depend on the Matsubara frequencies ujn = 
2TiT{n + 1/2), the position r, and the unit vector n on 
the Fermi surface, 

{un + n • [V - iA{r)]} f{uJn, n, r) = A{r)g{ujn,n, r) , 
{ujn - n-[V + iA{r)]} /(w„, n, r) = A^(r)g(a;„, n, r) , 

(1) 

where the dagger denotes complex conjugation. The 
(normal) Green's function g(u;„,n,r) is related to / and 
/ via the constraint (suppressing all the arguments for 
brevity) : 
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//=1. 



(2) 



These equations are to be solved together with the self- 
consistent equation for the complex order parameter 
A(r) and the Maxwell equation relating the magnetic 
field to the (super)currcnt: 



A(r)log|^ + 2^r^ 



Mr) 



,n,r) 



0, 



(3) 



V X H + i^2TTTy^ f—3ng(ujn,n, 



r) = 0, (4) 



with H = 'V X A. 

In Eqs. (HJ-lHl) we used as the unit of length the zero- 
temperature BCS coherence length 



2An 



(5) 



where vp is the Fermi velocity and Aq is the zero- 
temperature zero-field order parameter, which gives the 



energy unit. The vector potential A is rescaled by 
(/)o/27rfo and the magnetic field H by (/)o/27r^g, with 
00 = TTc/e as the flux quantum. These choices of units 
render all quantities dimensionless; for example, the BCS 
critical temperature is Tc — c'^^/tt ~ 0.567, where je is 
Euler's constant. Finally, the dimensionless parameter 
kq, the only independent parameter remaining after the 
units are chosen, is defined in analogy with the GL pa- 
rameter K as 



Ao 

where the zero-temperature penetration depth is 

'27reo^ ' 



1 

-^0 



(6) 



(7) 



with v as the density of states at the Fermi energy. 

In writing Eq. ([T]) higher-order terms in the mag- 
netic field are neglected which give rise to diamagnetic 
effects. This is a good approximation if 



(8) 



where ujc = eH/mc is the cyclotron frequency. This con- 
dition can be rewritten as: 



Af 1 H 
l^r^i?c(0) 



T 



(9) 



In low Tc, strong type II superconductors, the first two 
factors on the left-hand side are both small parameters. 
In what follows we consider magnetic fields H smaller 
than the zero temperature critical field ^^c(0); there- 
fore, this approximation is justified down to very low 
temperatures,^'* 



III. THERMODYNAMIC POTENTIAL 

Equations ([T]), ([3]), and (|4]) are the Euler-Lagrange 
equations obtained by varying the following functionaU^ 
with respect to / and /, At, and A, respectively: 



n^i^l d'r\ ^ {H{r) HaY + |A(r)|^ log ^ 



+ (dn) 



|A(r)|^ 



-At(r)/-/A(r) 



/ 



2tj„(5 - 1) - .gn • ( V log 4 - 2iA{r) 



(10) 



with g as implicitly defined by Eq. ([2]), Ha as the applied 
field which we assume uniform, and 
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(11) 
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The functional in Eq. ((T0|) is not the thermodynamic 
potential; however, for any given A(r) and A{r) Eq. pU]) 
gives the difference between the potentials in the su- 
perconducting and normal states once the solutions to 
Eq. Il]) for / and / are substituted into it. 

As a first step in solving Eq. we recall that the 
order parameter can be assumed as real; more precisely, 
the (gradient of the) phase of the order parameter can be 
collected together with the vector potential into a gauge 
invariant quantity. All other quantities become gauge 
invariant as well, and the new vector potential is propor- 
tional to the supercurrent velocity. This means that in 
the Meissner state A must vanish deep into the supercon- 
ductor; the same holds for its component perpendicular 
to the surface, as no current leaves the superconductor. 
Moreover, H — Ha at the surface. These arguments de- 
termine the boundary conditions for A; further boundary 
conditions are discussed at the end of this section. 

With a real order parameter, it is convenient to intro- 
duce the sum and difference of the anomalous Green's 
functions, 



s = f + f, d = f~f. 
Then, Eq. ([2]) can be rewritten as 



(12) 



(13) 



Taking the difference between the equations in Eq. ([T]) 
and solving for d in terms of s, we obtain 



d=- 



n ■ Vi 



where we introduced the short hand notation 
r2„ = Un — in ■ A . 



(14) 



(15) 



Our discussion so far is valid for any kq, but to find 
an explicit expression for the thermodynamic potential 
as a functional of A(r) and A(r), we look for an approx- 
imate solution to Eq. ^ for 3> 1 and at arbitrary 
temperature. To find a suitable approximate expression 
for s, we note that a rescaling r kqT only affects the 
gradient terms in Eq. ([T]), so that an expansion in the 
small parameter 1/ko is equivalent to a gradient expan- 
sion. To zero order we neglect the gradient terms, drop 
d in Eq. (jl3p , and using the sum of Eq. ([T]) arrive at 



„(o) 



2A 



(16) 



which for A = correctly reduce to the standard result*^ 
for a bulk superconductor in the absence of magnetic 
field. Hereinafter, due to the above-mentioned rescaling, 
Ao is the unit of length. 



To calculate the next order in the expansion, we define 
[see Eq. ^] 



1 n-Vs(") 



Kq fl„ 

From Eq. we obtain 

9 ^g^''^+9^'\ 



-,(2) 



,(0) 



.(2) 



(17) 



(18) 



^ 8g(0) 45(0)' 

where s'^^^ is the next non-trivial order in the expansion 
for s (i.e., s ~ s^'^^ -|- s^^-*), which is again found from the 
sum of Eq. H]), 



A(n-Vs(o))' 1 {n-Vfs 



(0) 



+ 



4 $72 Sn Sf-^ 

(n • V)i n ■ V.s(o) 



with 



(19) 



(20) 



It turns out that this expression, however, is not needed 
to obtain the thermodynamic potential, as its contribu- 
tions to it cancel out. This can be checked by substituting 
Eqs. (fT2|) - p3)l into Eq. (flO)) and using the approximate 
expressions in Eqs. P^ - (fT5|) . After an integration by 
parts (and dropping the resulting surface term), to low- 
est nontrivial order in the gradient expansion, the ther- 
modynamic potential as a functional of A and A is 



n 



(dn) 



(V X A-Haf + log 



A2 



2 



(21) 



_^yn|TA2 



n 



This expression is one of our main results and is the 
starting point to study the metastability of the Meiss- 
ner state; see Sec. IIVI It can be considered as an ex- 
tension of the GL approach, and as a check, we show in 
Sec, nil A1 that Eq. ((2T|) reduces to the known GL poten- 
tial in the appropriate limit. It is interesting to compare 
the present result with other extensions of the GL the- 
ory in the literature )^^'^°i^^'^^i^^ where the expansion is 
performed with respect to the covariant derivative [i.e., 
the gauge-invariant operator V — 2ieA^]. In the present 
notation, this amounts to supplementing the already per- 
formed gradient expansion with an expansion over A; at 
lowest nontrivial order the published results are recov- 
ered. As the field increases, however, A increases as well 
and this additional expansion is not reliable. In the lo- 
cal limit considered here the order parameter amplitude 
spatial profile is determined primarily by the depairing 
effect of the supercurrent, which is taken into account ex- 
actly, while the additional gradient terms are suppressed 
by the small parameter I/k^^ 
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A. Ginzburg-Landau limit 



The GL parameter kgl is proportional to kqi 



As remarked in the Sec. [H the GL approach is vahd 
near the critical temperature. More generally, near a 
second-order phase transition the order parameter is 
small and an expansion of the thermodynamic poten- 
tial in powers of the small parameter A/27rT becomes 
viable. To perform this expansion in the present case, 
we introduce the rescaled vector potential 



2 A 



3 A(r) ■ 



(22) 



Here A(T) is the value of the order parameter at tem- 
perature T in zero magnetic field, which by definition 
satisfies the equation 



1 



1 



(23) 



obtained by minimizing Eq. (|2ip with A — Ha — 0. 
The expansion of the above equation near Tc leads to the 
well-known^^ approximate expression 



MT) 
2-kT 



1 



T 



where 



(n + 1/2)3 



-^"( 2) = ^C(3). 



(24) 



(25) 



The rescaling in Eq. (j22p is equivalent to normalizing 
the field with respect to the temperature-dependent ther- 
modynamic critical field rather than its zero temperature 
value; the scaling is to be applied to the external field Ha 
as well. Similarly, we define the normalized order param- 
eter 



V'(r) 



A(r) 
A(T) 



(26) 



and introduce the temperature-dependent penetration 
depth 



A(T) 



Aq 2'kT Aq 



1 



(27) 



as the unit of length. It should be stressed that all the 
temperature dependencies introduced in this subsection 
are valid onl y in the vicinity of the critical temperature; 
see also Ref. [20- 

Substituting the definitions in Eqs. and ([25)) into 
Eq. (|2T|) . expressing lengths via A(r), and expanding in 
powers of A{T)/2ttT <C 1, the lowest-order term is 



a 



GL 



(27rT)2 



d^7 



I IV X A- H, 



(28) 



KGL = a/ ^ '«o ~ 0.42ko . 



It can be written as the ratio 



A(T) 



(29) 



(30) 



between the temperature-dependent penetration depth 
[Eq. (P7)) ] and coherence length 



C(T) 



2 2ttT 

3 A(r) 



^0. 



(31) 



Inside the curly brackets in Eq. ((28)) one can recognize 
the GL free energy in its dimensionless form; see, e.g., 
Refs. i and [^. The last term, in particular, represents 
the energy associated with the spatial variation of the 
amplitude of the order parameter. Due to this term the 
appropriate boundary condition at the surface is ip' — 0, 
where the prime indicates the derivative along the nor- 
mal to the surface. The similar term in Eq. (|2T)) has a 
more complicated dependence on both the order param- 
eter and the vector potential; for this reason the issue of 
the boundary condition for A demands further investiga- 
tion beyond the scope of the present work; for example, 
the addition of surface terms [such as those discarded in 
obtaining Eq. i)^] may be necessary; see, e.g., Ref. I27I. 
In what follows, however, we will concentrate on the limit 
kq ^ 00; in this case the last term in Eq. ()2ip is ne- 
glected, the order parameter is determined by a "local" 
equation (rather than a differential one - see Sec. IIV() . 
and no boundary conditions are required for A(r). 



IV. STABILITY CONDITION 

In Sec. mil we have derived an approximate expression 
[Eq. (I2T)) ] for the thermodynamic potential valid at large 
kq. a standard procedure can be now applied to study 
the properties of this functional; for example, the self- 
consistent equation for A and the Maxwell equation are 
found by taking the variations with respect to A and A, 
respectively. Considering from now on only the lowest 
order contributions in I/kq [i.e., neglecting the last term 
on the right-hand side of Eq. (|2T)) ] and using Eq. ((23)) . 
we find 



(dn) 



1 



1 



0, (32) 



which gives, via definition (fT5)) . a local relation between 
order parameter and vector potential, and 



V X V X A - i J{dn) 3n 



0. 



(33) 
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Solutions to Eqs. ((32)) and ([33]) are (meta)stable only 
if they are a minimum of Q; i.e., if its second variation is 
positive. To investigate the stability, let us parametrize 
A and A as 



A = As+T] , A = As + a, 



(34) 



where and As satisfy Eqs. p2|) and ([33)1 . Expanding 
Eq. ((2T|) for small 77 and a, the second variation is: 



-2i 



{dn) 



Al 



(f^2+A2)3/2 



(f]2+A2)3/2 

ri{n ■ a) 



I (V X af 



(35) 



with = ujn — in ■ Ag. In the absence of magnetic 
field (i.e.. As ~ 0) the superconducting state is stable 
for any T < Tc] after integrating over n the last term in 
square brackets in Eq. ([55)1 vanishes, so that S^ft > for 
any fluctuation as long as A^ > 0. As the field increases, 
however, the sign of S^fl changes, by definition, when the 
superheating field is reached. Therefore, to find Hgh we 
look for nontrivial fluctuations i] ^ 0, a ^ such that 
S^n[As,As] = 0. 

In the geometry under consideration (i.e., supercon- 
ductor in the a; > half space and Ha \\ z), the solution 
to Eqs. (|32p and ([55)1 is parametrized as 

A, = A,(a;), As ^ [0, Ay{x),0] . (36) 

Then, as shown by Kramer— in the GL limit, the fluctu- 
ations can be taken in the following form: 



(37) 



77 = Tji^x, k) cos{ky) , 

ct = [ax{x, k) sin(fcy), ay{x, k) cos{ky), 0] . 

and minimizing with 



Substituting Eq. ^7}i into Eq, 
respect to rj and dx, we find 

G 

V = Tray , 



(38) 



The first term on the right-hand side of Eq. (|42)) gives 
always a positive contribution to the second variation, 
but as we now argue it can be neglected in the large kq 
limit. Clearly, the larger k is the smaller this contribution 
becomes; on the other hand, the second variation of the 
last term in Eq. (PT|) . which we have neglected, would 
schematically contribute a (positive) term proportional 
to fc^/^Q. Therefore, the optimal value is fc '--^ which 
is in agreement with the GL result of Ref. 0, and the 
two terms both give contributions of I/kq, which we 
neglect as kq — )■ +00. As a consequence, the superheating 
field is determined by the vanishing of the coefficient of 
the second term on the right-hand side of Eq. (|42)) . i.e.. 



0. 



(43) 



In Sec. |V] we use this condition together with Eqs. 
and ([33]) to calculate the superheating field. 



V. SUPERHEATING FIELD 

The stability analysis of the previous section gives a 
conceptually simple procedure to find the superheating 
field; we should first solve Eqs. ([32]) and ((33)) to find the 
profile [cf. Eq. ([36]) ] of the (meta)stable superconducting 
state for a given applied magnetic field, and then find 
Ha such that the condition for the instability threshold 
in Eq. (|43)) is satisfied. The task of solving the non-linear 
differential equation (|33p. however, makes this route dif- 
ficult in practice. An alternative approach is based on 
the observation that the combination 



Hi 



3 [dn) 



2 



(44) 

where Hg = 'V x As, is constant throughout the 
superconductor,?^ as can be checked using Eqs. (|32p and 
(I33p . Taking into account the boundary conditions, the 
form of the solution in Eq. ([55]) . and the expression 



" 3F, + fc2^' 
with prime denoting differentiation with respect to x, 



(39) 



/", , iflsAsUy 



and 



= (dn) 



(f]2 + A2)3/2 



(40) 



(41) 



Here no = 1 and we note the property |_Fi| < 1. With 
these definitions, we obtain for the second variation (up 
to a numerical prefactor) 



S^n cx 



dx 



Fx 



3Fx 



Fo 



(42) 



H^{T)^67rTY^ 2(Va;2-hA(T)2 

n L 

_ A(r)2 

for the critical field, we find (cf. Ref. Isl) 



OJr 



(45) 



rr2 tt2 

ti„ — 



E 



Im ilo 



A2 



(46) 

where Aq — Ay{0), Qq = Un — iAo, and Aso = As(0). 
The above equation relates the applied field to the values 
of the order parameter and vector potential at the sur- 
face. At the superheating field, these two quantities can 
be found by solving the local Eqs. ([5^ and (H5|) . This 
can be done analytically in the limiting cases T — > Tc 
and T ^ 0, as we now show. 
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A. Limiting cases 

In the GL (T Tc) limit Eqs. §^ and ga]) reduce 
to, respectively, 



(47) 



4 

- 4^ A-^ 



Solving these equations we find A^o 
Defining 




H 



Hsh 



(49) 



and substituting the result into the GL limit of Eq. ((46|) 



= 1- 



A 



sO 



A{Ty 



(50) 



we arrive at H = V5/3 ~ 0.745, which IS m agree- 
ment with Refs. and |3- This is not surprising, since 
we showed in Sec. IIII Al the reduction of our thermody- 
namic potential [Eq. to the GL one [Eq. (gH])] hi 
this limit. 

In the opposite limit T — s- and using the notation 
A = A^oMo, Eqs. ^ and ^ become^ 



log (A 



sO) 



log [Ao(l + Vl^)] - - 

(A,o <Aq< e/2) 



(51) 



1 r 



1 - - A2 i 1 - - (1 + 2A2) 



while Eq. pS)) can be written as 



(52) 



= 0, 



i/^ = 1 - -Ai 



- - - (1 - A") 
2 2^ ^ 



1-A 



2n3/2 



(53) 



Substitution of the solution to Eqs. ((5T|) and ((52)) into 
the latter expression finally gives H ~ 0.840, as found in 
Ref. i. 



B. Temperature dependence 

Having verified that the known limiting results are re- 
produced with our approach, we now consider the tem- 
perature dependence of H. Near the two limiting tem- 
peratures we can in principle calculate temperature de- 
pendent corrections by further expanding over the small 
parameters A/2ttT close to Tc and T/Aq as T ^ 0. In- 
stead, to obtain the behavior at arbitrary temperature. 



0.84- 

0.82 - 

^ 0.80- 

0.78 - 

0.76- 

0.74 — 
0.0 



0.2 



0.4 0.6 

TIT, 



0.8 



1.0 



FIG. 1: Temperature dependence of the ratio Hgi^/Hc. Note 
the nonmonotonic behavior with a maximum at low temper- 
ature T fa 0.06Tc. 



we resort to a numerical approach. Following the same 
strategy as in Sec. IV Al we first solve numerically the sys- 
tem composed of Eqs. ([5^ and ([33]) to find and Ago] 
then, we substitute the result into Eq. ([46]) . In this way, 
we obtain the curve presented in Fig. [T] Interestingly, 
we find a nonmonotonic dependence of H on tempera- 
ture, with a maximum of H'^^'' « 0.845 at T « 0.06Tc. 
A nonmonotonic behavior is also found in Hsh{T) shown 
in Fig. [2] Taking into account the decrease in Hc{T) 
with increasing temperature, Hsh_{T) acquires its maxi- 
mum value m 0.843iJc(0) at the lower temperature 
T « 0.04rc, see the inset of Fig.E] Since i?2(0) = ^ttvAI 
depends only on material properties, this implies that 
in the London limit there is an optimal temperature at 
which the superheating field is the highest possible for a 
given material. 




FIG. 2: Temperature dependence of normalized by the 
zero-temperature critical field. The nonmonotonic behavior 
with a maximum at T « 0.04Tc is evident in the inset, which 
zooms in on the low-temperature region. 
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VI. IMPLICATIONS FOR ACCELERATOR 
DESIGN 

The rf cavities used in modern particle accelerators, 
both for high-energy physics and as x-ray sources, are 
made of superconducting niobium. The best Nb cavities 
are operated in the metastable region well above Hci- 
The superheating field Hgh provides an upper bound for 
the maximum particle acceleration that a given cavity 
can produce, and the operating point for the best cav- 
ities is approaching the theoretical limit provided by the 
GL theory, when the latter is extrapolated to the oper- 
ating condition (T 0-2Tc) where it is not valid. Nio- 
bium is not a high-K material, and our theory cannot 
be directly applied to it, but Fig[T] would indicate that 
at T/Tc ^ 0.2 the true superheating field would be 11% 
higher than the GL estimate for a high-K material. A 
change in the theoretical upper bound of this magnitude 
would have significant implications for future attempts 
to improve the material processing of existing Nb-based 
cavities. In principle, a numerical solution of the linear 
stability problem for the Eilenberger equations should be 
possible (albeit challenging) for all values of k, including 
K ~ 1 - a calculation of direct significance to current 
technological applications. 

There are several other superconducting materials 
which appear potentially promising as eventual replace- 
ments for Nb in future accelerator applications, all of 
which have significantly higher k and hence are poten- 
tially better described by our London limit calculation. 
For example, if run at the current operating tempera- 
ture of 2 K, cavities made of NbaSn or MgB2 would 
be near the peak of Hsh/Hc in Fig. [U and hence our 
calculation would suggest a peak field 13% higher than 
that provided by the GL theory. Using a current design 
for the superconducting cavity, our result for Tfgh sug- 
gests a theoretical upper bound for the accelerating field 
of 200 MV/m, a factor of 4 larger than the operating 
fields of the best Nb-based cavities. Material difficulties 
have so far kept high-temperature copper-oxygen-based 
superconductors from being useful in these applications, 
but new high-Tc materials, e.g., iron pnictides,— may 
provide more forgiving material properties. A quantita- 
tive estimate of the superheating field in these materi- 
als, however, may demand calculations that incorporate 
effects that go beyond the semiclassical analysis of the 
present work. For example, a complete description of 
superconductivity in MgB2 requires an Eliashberg-type 
calculation,™ with material properties extracted from a 
density-functional electronic structure calculatioui^ 

We point out that our result is in sharp contrast 
with the commonly-used heuristic Hsh ~ Hd n of Yogi 
et alM- This heuristic, termed as the "line nucleation 
model" , is not a linear stability calculation, but an en- 
ergy balance argument that gives a nonsensical estimate 
Hs\x < Hci for large k. The formula's success in de- 



scribing experiments^ - suggests that there may be nucle- 
ation mechanisms (perhaps disorder-mediated) that be- 
come more diflticult to control in high-K materials, but 
it should be viewed as an experimental extrapolation, 
rather than a theoretical bound, in guiding the explo- 
ration of new materials. 



VII. SUMMARY AND OPEN PROBLEMS 

In this paper we have revisited the problem of evalu- 
ating the superheating field for type II superconductors, 
in particular with regard to its dependence on tempera- 
ture. To extend previous calculations->^>^i^ based on the 
Ginzburg-Landau theory, which is restricted to temper- 
atures close to the critical one, we have employed the 
semiclassical approach in order to derive an approximate 
expression for the thermodynamic potential [Eq. (|2ip ] 
valid at large values of the Ginzburg-Landau parameter 
K ^ 1. From this expression we have calculated, in the 
limit K +00, the temperature dependence of the ra- 
tio between the superheating and critical fields which is 
presented in Fig. [TJ The relevance of our results to ap- 
plications in particle accelerators is discussed in Sec. IVII 

Many natural extensions of this work come to mind, 
beyond and in connection with those already mentioned 
in Sec. IVII For example, it would be interesting to 
study in more detail the critical variations, as done in 
Ref. y in the GL limit, and the finite k corrections. In 
our calculations we have considered the simplest possi- 
ble case, namely a clean superconductor with a spheri- 
cal Fermi surface. However, the semiclassical theory can 
easily accommodate anisotropics in the Fermi surface. 
The effect of bulk impurities can also be incorporated in 
the formalism. In superconducting cavities in the pres- 
ence of rf fields various mechanisms for the breakdown 
of superconductivity are associated with characteristics 
of the surface, e.g., the presence of surface impurities 
or steps caused by grain boundaries. ^« Moreover surface 
properties, namely specular vs. diffuse reflection, are 
known to affect the electromagnetic response of impure 
superconductorsj^ Therefore, it would be important to 
explore theoretically the impact of surface imperfections 
on the superheating fieldi^ 
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